Preprint typeset in JHEP style - PAPER VERSION 



UAB-FT-669 



Non-perturbative asymptotic improvement of 
perturbation theory and Mellin-Barnes representation 



S. Friot 

Institut de Physique Nucleaire d'Orsay 

Universite Paris-Sud 11, 91405 Orsay Cedex, France 

E-mail: samuel . f riotOipno . in2p3 . f r 

D. Greynat 

Institut de Fisica Altes Energies 

Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain 
E-mail: greynat@ifae.es 
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we show how to obtain hyper asymptotic expansions from the (divergent) formal power 
series which follow from the perturbative evaluation of arbitrary "A-point" functions for 
the simple case of zero-dimensional cfr field theory. This hyperasymptotic improvement 
appears from an iterative procedure, based on inverse factorial expansions, and gives birth 
to interwoven non-perturbative partial sums whose coefficients are related to the perturba- 
tive ones by an interesting resurgence phenomenon. It is a non-perturbative improvement 
in the sense that, for some optimal truncations of the partial sums, the remainder at a 
given hyperasymptotic level is exponentially suppressed compared to the remainder at the 
preceding hyperasymptotic level. The Mellin-Barnes representation allows our results to 
be automatically valid for a wide range of the phase of the complex coupling constant, 
including values for which the perturbative expansion is not Borel summable. A numerical 
analysis is performed to emphasize the improved accuracy that this method allows to reach 
with regard to the usual perturbative approach, and the importance of hyperasymptotic 
optimal truncation schemes. 
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1. Introduction 

The divergent behavior of a (divergent) asymptotic expansion does not at all detract from 
its computational utility. This statement is corroborated by the fact that, in what concerns 
its first few partial sums, a divergent asymptotic expansion of a given quantity "converges" 
in general much faster to the exact result than a convergent series representation of the same 
quantity does. In the case of the Standard Model quantum field theories, we may therefore 
even say that, regarding the extreme difficulty to go beyond the first few perturbative orders 
when computing observables in QCD or in the electroweak theory, it is an advantage, for 
phenomenology, to deal with a formalism that leads to presumably 1 asymptotic power 
series, diverging for all values of the coupling constants, rather than convergent ones. 

x It is not yet possible to prove that the formal perturbative expansions in the Standard Model quantum 
field theories are asymptotic expansions, but the fact that the sum of the first few perturbative terms is in 
general in very good agreement with "exact" experimental results gives a piece of evidence that perturbation 
theory is asymptotic to "something". The question is: to what? [1] 
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However, one has of course to keep in mind that when dealing with divergent asymp- 
totic perturbative power expansions, there always remains a finite limit of precision beyond 
which the usual asymptotic theory cannot go, even when the objects that one wants to 
compute are well-defined 2 . All ways to break open this precision limit are welcome. In 
the beginning of the 1990's, new asymptotic objects, which have in general a larger re- 
gion of validity (a larger domain of definition in the complex expansion parameter) and a 
greater accuracy than conventional asymptotic expansions, appeared in the mathematical 
litterature [2]. With them, a new asymptotic theory emerged: exponential asymptotics (or 
hyperasymptotics 3 ). These asymptotic objects (hyperasymptotic expansions) are very in- 
teresting since they correspond, for some optimal truncation schemes to be defined later, to 
what we could call in physics a non-perturbative asymptotic improvement of a perturbative 
(asymptotic) power series. Our aim in this paper is to show how hyperasymptotic expan- 
sions appear naturally in the simplest example one may have in mind in "particle physics", 
namely zero-dimensional (j) 4 field theory. In particular, we show that in the course of the 
study of the iV-point functions of this theory one may obtain hyperasymptotic expansions 
directly from the formal (divergent) expansions which follow from their perturbative eval- 
uation and without further information (i.e. without using the fact that by their integral 
representation, for instance, we have a rigorous definition for these objects). 

The paper is organized as follows. In the introductive section 2.1 where basic facts 
are recalled, the perturbative approach is detailed and, to fix ideas, numerical results are 
given for a particular value of the coupling constant A. The main part of the paper is 
section 2.2. In section 2.2.1 we present the formal approach which allows to rewrite the 
perturbative results in terms of hyperasymptotic expansions (at first level in the hyper- 
asymptotic process) . The calculations are based on so-called inverse factorial expansions of 
the ratios of Euler Gamma functions which constitute the coefficients of the perturbative 
terms forming the tail of the perturbative series. This makes appear the Mellin-Barnes 
representation into the game. After term by term Borel resummation of the reexpanded 
tail (this strategy is inspired from [5]), a new expansion of the tail emerges which, added to 
the perturbative partial sum, form the hyperasymptotic expressions at first hyperasymp- 
totic level. At this level one may already notice a resurgence phenomenon that links the 
perturbative coefficients with those of the tail's new expansion, and which will also be 
observed at each higher hyperasymptotic level. Next, section 2.2.2 explains in detail how 
the hyperasymptotic theory of Mellin-Barnes integrals [4] leads to the proof that the hy- 
perasymptotic results obtained in section 2.2.1 are correct, and this for a wide range of 
the phase of the complex coupling constant (taking into account, this time, the integral 
representation of iV-point functions that was avoided until here), in particular on Stokes 
lines where the perturbative expansions are not Borel summable. In section 2.2.3, one is 

2 In the case of the (4D) Standard Model of particle physics, due to the absence of a definition of the 
theory, already the correct evaluation of the size of this precision limit is lacking (i.e. theoretical errors 
implied by truncations of perturbative expansions, OPE, etc. are not under control) although it is in 
principle of crucial importance in precision test of the Standard Model if one aims to find new physics 
effects. 

3 The first paper dealing with hyperasymptotics is [3], but here we mainly rely on [4] where the theory 
is developed for Mellin-Barnes integrals. 
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concerned with optimal truncation schemes of the hyper asymptotic expansions at zeroth 
and first hyperasymptotic levels. Their link with the non-perturbative interpretation of 
our results are underlined. Then follows section 2.2.4 where higher order hyperasymptotic 
levels are obtained. One ends the main body of the paper in section 2.2.5 by performing a 
numerical analysis which allows to compare the hyperasymptotic expansions for different 
optimal truncation schemes with the perturbative results. Two appendices give the proofs 
of some results quoted in the text. 

2. Exponential asymptotics in zero-dimensional Euclidean </> 4 theory 

The 4-dimensional Euclidean (ft 4 action is given by 



S = d A x 



U^{x)d»(ft{x) + \m 2 <f{x) + j,<ft\x) 



(2.1) 



Going to the O-dimensional theory, which consists to reduce space-time to just one point, 
makes that the x-dependency disappears (4>(x) becomes a simple real variable (ft, the deriva- 
tive term and the overall integral disappear), so that one gets 

S(<ft) = \m 2 <? + ^<ft\ (2.2) 

The corresponding 4-dimensional generating functional is therefore also reduced, in the 
O-dimensional case, to a usual integral of the form 

d<ft e -§-V-^ 4 +^, (2.3) 

where j is the external source and M a normalization factor. 

In the following, objects under study are the "iV-point" functions, defined as 



where 



deft e -2™V- a 4 (2.5) 

-oo 

is the O-dimensional version of the generating functional of vacuum to vacuum transitions. 
Notice that 

G (2p+1) = 0. (2.6) 

Eqs. (2.4) (and (2.5)) are defined for Re A > and we want to study the small A 
case. We will see that, thanks to the Mellin-Barnes representation, the results that will 
be obtained in this paper are in fact valid for a much wider range of complex values of A 
than just the right half complex A-plane. In particular, our results are valid on Stokes rays 
where the perturbative expansion is not Borel summable (in our case | arg A| = ir). 

Since the theoretical results of interest that we found may already be observed in the 
study of Z(0) with N = we think that it is more pedagogical to present detailed 
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derivations of our analysis on the simple example of Z(0) rather than on arbitrary iV-point 
functions. This avoids a dependance on the iV parameter in the calculations, which would 
complicate the expressions in the text without really giving new results compared to Z(0) 
(apart from the form of the resurgence phenomenon, see section 2.2.4). However, the final 
expressions of the calculations will also be given in the more general Appoint functions 
case and some subtle changes that have to be done in the computation, for our theoretical 
strategy to be valid in this more general case, will also be explained when necessary. 

2.1 Pert ur bat ive approach 

In this introductive section, we recall basic facts about the link between the perturbative 
expansion of O-dimensional Euclidean (f) 4 theory and the counting of Feynman diagrams in 
the corresponding 4-dimensional theory. We do this on the example of Z(0) and, to fix 
ideas before the non-perturbative asymptotic analysis, we also give perturbative numerical 
predictions for a particular value of the coupling constant. 

2.1.1 Perturbative expansion 

So our first object of study is Z(0) with 4 m = 1 and M = 

Z(0) = / dtp e-W-v**. (2.7) 

V27T J-oo 

In this very simple case, one may get a closed expression for (2.7) in terms of special 
functions or even keep the integral representation and perform the whole (non-perturbative) 
asymptotic analysis in a fully rigorous way. This will be done as a check, but what we 
precisely want to show in this paper is a method for performing the analysis in the con- 
verse way. Indeed, the question we want to answer is: how can we get a non-perturbative 
improved expansion from a divergent perturbative expansion when only the latter is avail- 
able? 

Let us compute the perturbative expansion of (2.7) for small A with Re A > 0. Re- 
placing e ^ by its series representation and performing term by term integration, one 
finds 

For further purpose, we define u& so that Z(0) ~ ^2t=o u k an< ^ the perturbative partial 

sum^r = £Lo^- 

The first few terms of (2.8) are given by 

Z (0) = 1 - I A + ^A 2 - -^A 3 + ^A 4 + 0(A 5 ). (2.9) 
w 8 384 3072 98304 v ; y J 



4 We could keep the mass as a free parameter but putting it to unity, as well as choosing the normalization 
constant of Z(0) to be (2n)~2 , makes the counting of the 4-dimensional Feynman diagrams, that the 0- 
dimensional theory allows, more transparent right from the beginning. 



- 4 - 



The expansion in the r.h.s of (2.8) is divergent for any value of A, as can be seen from 
the fact that 



Uk+l 



Uk 



(_i)*+i r(2(fc+i)+|) fc+1 

(fc+l)! (3!) fc +! A 



(-i) fc r(2fc+i) 

~S1 (3!)^ A 



1 (2fc + §)(2A; + ±; 



3! 



->• +oo, 



(2.10) 



and it is easy to prove that it is an asymptotic expansion of Z(0) as we will see later, but 
for the moment we only suppose that it is so 5 (we do as in the case of Standard Model 
gauge theories where one assumes that perturbative expansions are asymptotic to objects 
whose definition is still lacking today). 

The perturbative result for an arbitrary iV-point function (with N = 2p) is similar 



since 



+oo 



_1^2_ A ,4 



■1)' 



fc=0 



k\ 



(2.11) 



2.1.2 Feynman diagrams counting 

It is clear that the coefficients in the expansion (2.8) are related to the counting of Feynman 
diagrams in the 4-dimensional theory (taking into account their symmetry factors), since 
at first order in the A-expansion we have, for 4-dimensional Z(0), the contribution 



(2.12) 





and, at second order, 



(-1) 2 6x6x2 ,„ 



2! 



(4!) 




(-1) 2 3 x 3 2 
+ 2! (4!) 2 



+ 



(-1) 2 4! 
2! (4!) 2 




24 
384 




(2.13) 



where, as an example of how we get the coefficients in the l.h.s of (2.13), the 4! in the 
numerator of the coefficient in front of the last diagram is the number of different Wick 
contractions of the fields that lead to this topology (see Figure 1). 

We see that by summing the coefficients of each different topology (and this defines 
our counting of Feynman diagrams), one obtains from (2.12), (2.13) and higher orders the 



3 This is why we wrote an asymptotic equality in (2.1 
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Figure 1: The 4! different Wick contractions leading to the topology of last diagram in (2.13). 



result written in (2.8), see also (2.9), apart from the first term that does not exist in the 
4-dimensional theory, since there are no tree vacuum diagrams. 

Zero-dimensional field theories therefore have a practical interest for 4-dimensional par- 
ticle phenomenology by the fact that they allow a partial but important check of Feynman 
diagrams coefficients appearing in perturbative calculations. 



2.1.3 Numerical analysis 

Let us now perform a bit of numerical analysis to see the efficiency and predictive power 
of the perturbation theory developed in the preceding subsections. In the following, we 
choose A = \. 

Imagine that we do not know the value of (2.7) for A = | so that the only information 
that we have for computing Z(0) is the divergent perturbative expansion in (2.8) (remember 
that we made the hypothesis that it is an asymptotic expansion of Z(0)). Let us see what is 
the best perturbative prediction that can be obtained from (2.8). In Table 1, we computed 
the first twenty truncated sums 

iji G {1, 20}) of (2.8) for A = \ (and the values of 
the general term of the perturbative series for the same values of n). 

It is readily seen that the partial sums rapidly converge to a value around 0.965 from 
which they finally begin to diverge (see also Figure 2). 

In fact, since (2.8) is a divergent alternating asymptotic series, simple general inter- 
pretative considerations for its sum lead to the fact that, if we define the remainder R n so 
that 

ST = S^\ + R n , (2.14) 
one is driven to (see Appendix A for details) 



\Rn\ < 



\Ur. 



(2.15) 



and 



\Rn\ < \u n -\\ 



(2.16) 
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k 


Uk 


n 


OPcrt v~^v n — 1 




k 


Uk 


n 


C*Pcrt v~^v 71 — 1 

a n-l ~ 2^k=0 u k 





1.00000000 


1 


1.00000000 




10 


0.02354141 


11 


0.98145758 


1 


-0.04166667 


2 


0.95833333 




11 


-0.05240343 


12 


0.92905416 


2 


0.01012732 


3 


0.96846065 




12 


0.12827922 


13 


1.05733338 


3 


-0.00464169 


4 


0.96381896 




13 


-0.34248906 


14 


0.71484431 


4 


0.00314281 


5 


0.96696177 




14 


0.99043216 


15 


1.70527648 


5 


-0.00281980 


6 


0.96414197 




15 


-3.08409571 


16 


-1.37881923 


6 


0.00315269 


7 


0.96729466 




16 


10.2883505 


17 


8.90953129 


7 


-0.00422235 


8 


0.96307231 




17 


-36.6060184 


18 


-27.6964871 


8 


0.00659010 


9 


0.96966241 




18 


138.374139 


19 


110.677652 


9 


-0.01174624 


10 


0.95791617 




19 


-553.800008 


20 


-443.122356 



Table 1: Numerical values of the perturbative general term and partial sums of (2.8), for A = |, 
with an 8 decimal places precision. 
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1 ■ 1 

0.98 ■ 0.98 

0.96556048 

0.96 . . 0.96 

0.94 0.94 

0.92 0.92 

1 3 5 7 9 11 13 



Figure 2: S^-l (n S {1, 12}), for A = |, compared to the exact result (2.20) (horizontal line). 

We may conclude from these inequalities that the best truncation of the series (obtained 
by minimizing the remainder R n ) is theoretically obtained by truncating before or after 
the smallest term in magnitude. 

Therefore, there are two "best predictions" from the perturbative expansion (2.8) with 
A = ^, which only differ by their central value (see Table 1): 

Z(0)|flV = 0.96696177 ± 0.00281979 (2.17) 

A— 3 

and 

Z(0)|^i 2 = 0.96414197 ± 0.00281979. (2.18) 

It is clear that since (2.8) is an alternate series, the exact value has to be between them, 
from what we conclude that 

Z(0)|£?i = 0.96555187 ± 0.00140990 (2.19) 
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so that the perturbative expansion (2.8) leads to an already very good prediction (from 
the precision level viewpoint), of the order of 0.15%. 

In fact, from (2.7), one gets, with an 8 decimal places precision, 

Z(0)| A= i = 0.96556048. (2.20) 

The central value in (2.19) is therefore very close to the exact value, and corresponds actu- 
ally to the standard Stieltjes approximative resummation formula for alternating divergent 
series, which reads [5] 

r,-l 

Z(0)~J>* + 2^' (2 ' 21) 

fc=0 

where u v is the term of least magnitude. 

The point is that one can approach the "right" value (2.20) much closer than what 
perturbation theory does, by a refined asymptotic analysis that we present in the next 
sections where we obtain non-perturbative asymptotic improvements of the perturbative 
expansions (2.8) and (2.11). At the first stage of this analysis, the improvement takes the 
form of so-called exponentially improved asymptotic expansions [2] but, the process being 
iterative, it will be possible to get more and more non-perturbative refinements, in terms 
of hyperasymptotic expansions of higher level. 

2.2 Exponential asymptotics: non-perturbative asymptotic improvement 

To go beyond perturbation theory, one could perform a Borel resummation of (2.8) and, in 
this simple example, it works: one may reconstruct Z(0) from its perturbative expansion 
by a Borel resummation, for Re A > 0. 

On the one hand, as we will see, the method we present here also allows to reconstruct 
Z(0) from its perturbative expansion but, on the other hand, it gives another look at Z(0): 
as a non-perturbative asymptotic method that reveals interesting effects that would be 
hidden in a usual Borel resummation, like a resurgence phenomenon. 

2.2.1 Interpretation of the divergent perturbative expansion 

Numerically, we saw that the first few perturbative terms already do a very good work for 
the description of Z(0), but that the divergent character of the perturbative expansion is 
inavoidable if one includes more and more terms. In order to solve this problem, we have 
to give a meaning to the tail of the perturbative expansion. With this in mind, we divide 
the perturbative expansion into two parts, following [5]: 

j ~° ~^h,~- — w~ x+ ^hr* — <^ (2 - 22) 

The first part is the perturbative contribution that one wants to keep and the second, the 
tail of the divergent series. 

We are now going to perform formal manipulations that will be justified a posteriori, 
in the next section. 
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First, it is convenient to use the duplication formula 

1 



T(2z) = ^=2 2z - 1 T(z)T[z + ^) (2.23) 



to rewrite the tail as 

i ^ (-i) fc r(i + 2 fc) i ^ r(fc + i)r(fc + l) /-2A\ fc 

Our main tool is the so-called inverse factorial expansion which may be obtained from 
Barnes Lemma [6], [4] 



1 + , 4 / 1 4j = Y (-l) j AjT(k + — ds f(s)T(k - s) , (2.25) 

r(A; + l) } 3 K J> 2mJ c+m „ loo 1 > 



where 



r(i)r(|) 



, _ r( J + i)r(j + f) 

3 " 1W' ( } 



and the contour in the Mellin-Barnes integral on the r.h.s. is a straight line with c £ ] — 1, 0[ 
(for m = 0, c £ ] — |, 0[ but we always take m > 0). It is important to note that (2.25) is 
an exact equality 6 only for c + m < k (since in our case m and k are integers, and since 
min(fc) = n and c < 0, this is equivalent to m < n), see Appendix B. 

Inserting (2.25) in (2.24), exchanging the two sums and the sum and integral signs, 
one finds 



+ ^^L-,. ds/(s) £ r(t - s) (-) • (2 - 28) 

Now we perform Borel resummations using the definition [5] 

+ °° x -2X\ k „, , /-2A\ n . / 3 



k=n 



E r <* - i) ( V ) = r <" - « ( V ) A — ( n ) (2 - 29) 



where A^(x) is one of the so-called terminant functions [5], defined (when Re£ > — 1 and 
| argx| < 7r) 7 as 

1 f°° v e e" y 
6 One may in fact also prove that it is an asymptotic equality [6]. 

7 This terminant function can be extended to Ret < —1 via, for instance, its absolute convergent expan- 
sion [5]. 
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which can also be expressed as 



A e (x) = x l+1 e x T (-£, x) , (2.31) 
where we recall that F(a,x) is the incomplete gamma function defined, for | argx| < it, as 

/>oo 

T(a,x) = / dyy^e-y. (2.32) 

J X 

Notice that the terminant is nothing but the Borel integral of a general UV renormalon 
pole (see for instance [7]). It may also be seen as a Mellin transform. 
At the end, one obtains (for m < n) 

i ^ (_gr(|+2fc) 

k\ (3!) fc 

k=n 

= g(-i)^r<n - i) (A) ~ J r (_„ + , + J.) 

(-*\ \ n s 1 rc+rn+ico / o \ -s / o \ 

+ ht^L-^ d - (s) «»> r( - - s)r (- + « + 2a) • < 2 - 33) 

We therefore conclude that, for m < n and | argA| < it, 

7(a)- 1 y(-i) fc r(i + 2fc) fc 

m— 1 / „ \ -•)' / g 



r[-n-M + i, 2A 



^-g(-W-,)(A)- J r (- 

C_1l" •> 1 /-c+m+ioo / o \ -s / o \ 

+ ht '"m L.» d ' (s) «»» r( " - s > r (- + 5 + ^ s) ■ < 2 - M) 

The tail of the divergent perturbative series (2.33) has been rewritten as a partial 
sum, supplemented by a remainder integral written as a Mellin-Barnes representation. We 
therefore converted an infinite sum into a finite sum plus a convergent integral and it will 
be proven in the section 2.2.2 that the formal expression (2.34) is exact for any n and 
m as long as n > m and constitutes what is called the first level of the hyperasymptotic 
expansion of Z(0). Better than that, we will also prove in section 2.2.3 that the r.h.s of 
(2.33) is, at an optimal value of n to be defined later, exponentially suppressed with respect 
to A, so that (2.33) gives in fact, for this optimal value of n, the expression of a purely 
non-perturbative quantity. 

Eq. (2.34) is valid in a wider sector (| argA| < tt) than the usual Borel resummation 
which is valid for | argA| < ^. Moreover, we would like the reader to notice the similarity 
between the perturbative partial sum (first line of (2.34)) and the partial sum in the second 
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line of (2.34). Indeed, choosing n = 5 and m = 5, we explicitly get from (2.34) 

. . A 35 , 2 385 . 25025 . 

Z(0) = 1 - - + A A 3 + A 4 

v ' 8 384 3072 98304 

I_ e A { r(5)r (_ 4 , JL) _ | A r(4)r (_ 3 , A) + » * r(3)r (-2, A 



3072 



V 7 V 2 V 98304 v ; V 2A / J 



***IZ ds (^) /(-)r(5- a )r(-4 + a ,-). (2.35) 



V27T 2m j c+5 _ ioo v 2 V V ' 2A , 

We see that the coefficients of the second partial sum are the same as the perturbative ones 
(up to Euler gamma functions that we wrote explicitly to emphasize the symmetry of the 
formula). We chose n = 5 because, as we saw in (2.17), it is one of the two best orders for 
truncating the perturbative series for A = |. Let us however underline that this interesting 
phenomenon is independent of the choice of n. Indeed, for example, an equivalent formula 
to (2.35) is 

Z(0) = 1 - - + —A 2 
v ; 8 384 



_3_ 

e2A 



^ { r < 3 > r (- 2 - JO - ^ F < 2 » r JO + SS* ^ (»• 2A 



1 t 1 fC+3+ioo / o \ s / o \ 

W/^ L_ i00 ds (») m F(3 " s) F (" 2 + ») • (i36) 



r-c+3+ioo 

TT_ / 

where one sees that although the terms in the second line of (2.36) are not the same as those 
in the second line of (2.35), their coefficients are still equal to those of the perturbative 
contributions. This is a so-called resurgence phenomenon [4], which will also manifest 
itself at higher order in the non-perturbative asymptotic improvement process. It can 
be understood, at our present level, from the fact that the perturbative series (using the 
duplication formula and noting that y/2ir = r(^)r(|)) may be rewritten as 

^ C» 7? £ £ W^ At ~ 2'- 1 ''- 4 * (JO " ■ (2 ' 37) 

v fc=0 v ; fe=0 v 7 

which has the same coefficient than those appearing in the second line of (2.33). 

An alternative expression for Z(0) (completely equivalent to (2.34)) is obtained by 
inserting in (2.30) the Mellin-Barnes representation 

d+ioo 



[ A (yy* 71 (238) 

J 2m \x) sin(vrt) ' k ' ' 



l + a 

X 

d—ioo 

where d = Ret G]0, 1[. Performing the y-integral one then finds 

A e (x) = —J-— [ — (-) -^—T{l + l-t), (2.39) 
v ' r(£+l) J 2iir \x J sin(vrt) v h y > 

d—ioo 
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with Re (I — t) > — 1, so that, at the end, 

n-l , „ N -k 

Z(0) = ^(-l) k A k - 

k=0 ^ 

- — / dsf(s)— dul^r) . ? T(u-s) (2.40) 

V^vr 2?vr y c+m _ ioo 2?vry e+n _ ioo V 2 -V sm(vrn) 

with n > m, c €] — 1,0[ and d g] — 1,0[. In the double Mellin-Barnes integral, e G 
] — min(l, n — m — c), 0[. 

For the iV-point function case (with N = 2p > 2), our formal strategy applied to (2.11) 
imposes an inverse factorial expansion whith constraints on n and m different from those 
implied by (2.25), because of the p dependance of the ratio of Gamma functions. 

Indeed, the inverse factorial that we need is 

r(fc + § + l)r(fc + f + f) 
r(fc + i) 

m ~ 1 i /-c+m+ioo 

= ^(-r/£f r(fc + p-i) + — / ds5p ( s )r(fc + P - s ) , (2.41) 

■_q ^ 7r J c+m—ioo 

where the contour of the Mellin-Barnes integral in the r.h.s. is a straight line with c a real 
number so that c + m > — j + f and c + m < n + p, and 

_ r(.-S + j)r(.-j + }) r(-.) 

*W- r(i-|)r(l-f) (2 ' 42) 

and 

b (p) _ r(j-§ + |)r(j-§ + |) 

J " ;'r(i-i)r(i-i) • (2 - 43) 

Contrary to (2.25), one cannot naively put m = in (2.41). Indeed, Barnes lemma has no 
fundamental strip in this case since p > 0. To be valid here, Barnes lemma needs a deformed 
contour of the Mellin-Barnes integral which separates the poles of T (s — | + |) T (s — | + |) 
from those of T(— s)r(fe + p — s), while going from — ioo to ioo. 
Taking into account these facts one finds, at the end, 

2* / 3 V a .<_r„i_. , / 3V J ' / 3 



+ <- 1 >"S(^) % * E(-i)^(»^-.')(|;)" 3 r(-„- P + J + i ,A) 

c+m+ioo 

, ,J 2 V3V ^ /■ ds ( 3 \~ s , , , N / 3 \ 



c+m—ioo 

(2.44) 
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2.2.2 Mellin-Barnes hyperasymptotic theory: first level 

As we shall see now, the result (2.34) can be entirely justified from the modern point of 
view of Mellin-Barnes hyperasymptotic theory [4] if one does not anymore ignore, contrary 
to what we did since the beginning, the integral representation of Z(0). Indeed, after 
having computed the asymptotic expansion of Z(0) from its Mellin-Barnes representation, 
it will be possible to obtain from the expansion of the Mellin-Barnes asymptotic remainder 
integral, the hyperasymptotic expansion of Z(0) at first hyperasymptotic level. As a by 
product, it will be easy to show that (2.40) is valid for a wider range of phase than the 
complex A-plane with a cut on the negative real axis. In particular it will be proven to be 
valid for this ray (the negative real axis) on which (2.8) is not Borel summable (and which 
is a Stokes ray). 

Let us begin the calculations. Using the Mellin-Barnes representation 

j rc+ioo / ^ 



»r^*u*; r(s) ' (2 - 45) 

which is valid in the semi-infinite fundamental strip defined by c = Res e]0,+oo[ (the 
right half-complex s-plane) and for | argA)| < §, one obtains, from (2.7), 

-i r-c+ioo / \ \ —s r+oo 

Z(0) = ^L/ S U) ns)J_^r»*-**. (2.46, 
The 4> integral can be computed and we find 

where c = Res G]0, \ [ so that the (f> integral reduced the semi-infinite fundamental strip to 
a finite one. Moreover our Mellin-Barnes representation increases the A-domain of validity 8 
of Z(0) to | argA| < 4^, therefore (2.47) consists in an analytic continuation of (2.7) on the 
Riemann surface of the logarithm (we recall that (2.7) is defined in the right half complex 
plane only, i.e. for | arg A| < §). 

Now that we have this simple Mellin-Barnes representation for Z(0), it is straightfor- 
ward to get its asymptotic expansion when |A| — > 0. 

Indeed, for a closed rectangular contour, Cauchy's theorem gives (with n a positive 
integer and 8 e]0, |[) 



rc+ioo r-c—n+8+ioo r-c—n+S—ioo re— too 

+ + + 

J c—ioo J c+ioo J c—n+6+ioo J c—n 



ds A- r(s)r «- 2S > 



H'K^pK J c—ioo J c+ioo J c—n+8+ioo Jc—n+8—ioo_ (3!) 

- 1 2:TYV- 1)fc r( * + 2fe > (2 48) 



8 The A-domain of convergence of a Mellin-Barnes integral of the type (2.47) is given by |argA| < 
f (N n — Nd), where N n is the number of Gamma functions in the numerator of the integrand and Nd 
the number of Gamma functions in the denominator, taking into account the multiplicity of each Gamma 
function, the multiplicity of a Gamma function being defined as the absolute value of the number multiplying 
the variable in the argument of this function [4] . For instance, in the case of (2.47), we have N n — 1+ 1 — 2| — 3 
and N d = 0. 
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which is equivalent to 
Z(0) 



i g (-i)fcr(j + 2fc) Afc | 1 



v fc=0 



(3\) k 



ds A~ s V , 9 \ 2 _ 9 S (2.49) 

<i— n— ioo 



(3!)- 



where d = c + 5 e]0, 1[. Indeed, the integrals with paths f^ + - 



c— n+S+ioo 



c+ioo 



and J 



c—n+5—ioo 



are 



zero due to the exponential decay of the modulus of the integrand when the imaginary part 
of s = a + ir goes to infinity (r — > ±oo) and c— n + 5<a<c. 

This can be proven by using Stirling formula, which says that, for r — > ±oo, 



|r(<7 + »T)| = 0(|r| CT -^e-fl r l) 
One therefore has, for A = |A|e* 9 , 

r( s )r(±- 2s ) 



(2.50) 



(3!)" 



r^±oo 



(|Tr*e-iM)o(|T|* 
O (|r|- ff -§ e -¥M e ^) 



r| )|(|A|« 



id\—a—iT 



(2.51) 



The exponential decay is then obtained for |0| = | argA| < 4^. 
From now on we define the remainder integral of (2.49) as 

1 



R, 



2in^ J c - n _ ioo (3!)- s 



so that in the notation of section 2.1.1, 

Z(0) = S p n -\ + R n . 



(2.52) 



(2.53) 



To show that the exact relation (2.49) (or (2.53)) is an asymptotic expansion when |A| — ► 
(valid, as we will see, for | arg A| < 4f), one has to prove that 



R„ = 0(X n+a ), 



(2.54) 



where a > —1. Using the fact that \ J f(z)dz\ < f \f(z)\ dz, and performing the change of 
variable s = d — n + it, we find 



|R "'-UJ dt€ 



T{d-n + it)T [ - - 2(d - n + it) 



(2.55) 



This proves (2.54) by definition of the O symbol with a = —d, since the integral in the 
r.h.s of (2.55) is independent of |A| and converges for \9\ = |argA| < \ (see (2.51)). 
Therefore, 



Z(0) 



-i)*r(i + 2fc) 



(3! 



(2.56) 



which agrees with (2.8). 

Notice that the Mellin-Barnes procedure developed in [8] to compute Feynman dia- 
grams (see also appendices of [9] for details, as well as [10] and [11]) would have given 
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the same result (2.56), but only for A positive. This is due to the fact that the converse 
mapping theorem given in [11] has not been established for the complex A case. 

Recall that one wants to prove eq. (2.34). One starts from (2.49). 

Eq. (2.34) will in fact emerge from an expansion of the remainder integral R„ defined 
in (2.52). 

For practical reasons, it will be easier to perform our analysis if we rewrite (2.52) by 
making the change of variable t = — s and by using the duplication formula (2.23). 
This gives 

r - - i^k £T dt (s) " r( -" r H) r H) • < 2 - 57) 

where now d = Ret g] — 1, 0[. 

The hyperasymptotic approach developed in [4] for Mellin-Barnes integrals lies on 
inverse factorial expansions of the integrands. To apply it on (2.57), we use the reflexion 
formula 

iiKM = r(4)r(1 -«> < 2 - 58) 

to get 

r = L_ r^dtj^Y - r(t + |)r(t + t) 

2my/2ir Jd+n-ioc V2Ay sin(Trf) F(t + 1) K ' J 

The quotient of Gamma functions can be expressed as the inverse factorial expansion (2.25). 
One therefore has 



2m Jd+n-ioc V2Ay Sin(TTt) V J > 

4 — 7-^ / dsf(s)T(u-s), (2.60) 

V2A; sin(Trn) 2mJ c+m _ iao h K ! 



r- , -I du\ — 

V27T 2Z7T Je+n^joo 

where d e] — 1,0[, e g] — 1,0[, c g] — 1,0[, | argAj < 4^ and the coefficients are given 
by (2.27). The constraint from the inverse factorial expansion is Reu > c + m which is 
equivalent to e > c + m — n. 

This ends here the first level in the (rather simple) hyperasymptotic process. 

Now, this last formula is exactly (2.33) if one replaces the incomplete gamma function 
in the second line of (2.33) by its Mellin-Barnes representation 

r ( — n + x + 1, — 



2A 



-1 



/ o \ -n+x a 1 fd+ioo / o \ -t 

— e"2A / rft _ Y(t + n-x), (2.61) 

V2A; 2i7r7 d _ ioo V2Ay sin(Trt) V J> K ) 



r(n - x) 



with ci = Ret G] — 1,0[, in the third line of (2.33) by the same formula but with d = 
Ret g] — min(l, n — m — c), 0[, if one performs the change of variables t' = t + n and if 
one exchanges integrals in the second term of the r.h.s of (2.33). Equivalently, starting 
from (2.40), exchanging the two integrals in the third line of (2.40), one also get (2.53) 
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with R n given as in (2.60). Another way to prove it is to simply recover (2.59) from (2.33). 
This may be done by putting m = in (2.33). This is in fact also equivalent to restart 
from (2.24) with the m = version of (2.25) which is nothing but Barnes Lemma, see 
Appendix B. We have therefore proven that the interpretation of the tail of the divergent 
perturbative expansion given by (2.33) is correct and that (2.34) constitutes the first level 
of the hyper asymptotic expansion of Z(0). 

An important remark is that, written in terms of Mellin-Barnes representations, our 
expansion (2.40) obtained initially for values of A in the complex plane with a cut on the 
negative real axis, is now automatically valid in the wider sector | argA| < ^ (see footnote 
8) where it also gives much more precise results than (2.56). 

Let us illustrate this fact by computing Z(0)\ x= _i. From (2.47), after performing the 
change of variable t = —i(s — c), one gets numerically, with an 8 decimal places precision, 

Z(0)\ x= _i = 1.05995021 - 0.007584721 (2.62) 

Now, from (2.56) which is valid for A = — | but not Borel summable, the best prediction 
is 9 

Z(0)| A= _i = 1.06098837 ± 0.00140990. (2.63) 

so that the purely perturbative approach of course completely misses the imaginary part. 

In fact, imaginary contributions appear from the second sum of (2.40). Indeed, trun- 
cating the perturbative series after the fifth term and including contributions of the first 
two terms of the second line of (2.40) we get 

Z(0)| A= _i = 1.05990083 - 0.007527941 (2.64) 

We see that even without taking into account the remainder integral (third line of (2.40)), 
whose contribution would make us fall on the exact result (2.62), we obtain a very good 
description of Z(0) also for A < 0, contrary to (2.56) wich is not Borel summable for this 
value of the coupling constant. 

2.2.3 Non-perturbative asymptotic improvement of perturbation theory: opti- 
mal truncation schemes 

In section 2.2.1, we formally gave a meaning to the tail of the divergent perturbative 
expansion of Z(0) and of an arbitrary Appoint function, in terms of a partial sum added 
to a Mellin-Barnes remainder integral. This allowed us to obtain eqs. (2.34), (2.40) and 
(2.44) which are in fact hyper asymptotic expansions at the first hyper asymptotic level. 
Our interpretation of the tail has been rigorously justified in section 2.2.2 and we saw in 
(2.64) that numerically the partial sum in the second line of (2.40) gives contributions that 
perturbation theory (2.56) does not see (in particular an imaginary part, see (2.63)), so 
that of course there is something beyond perturbation theory in (2.40). In this section we 
explain the non-perturbative nature of these results. The crucial point is that "by allowing 
the number of terms in an asymptotic expansion to depend on the asymptotic variable, it 

9 Eq. (2.63) is obtained in the same way as (2.19). 
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is possible to obtain an error term that is exponentially small as the asymptotic variable 
tends to its limit" [2]. In other words, there is an optimal value of the truncation index 
of the perturbative series for which the remainder of the expansion is a non-perturbative 
quantity and in our example this optimal value is obtained by minimizing the remainder 
with n scaling like -Xr. This will define optimal truncation schemes. Let us see this in more 
detail. 

The perturbative expansion of Z(0) leads us to (2.49) where the integral in the r.h.s. 
may also be written as (2.59): 



R, 



1 



d+n+ioo 



2ilT\f2ir Jd+n-ioo 

where d = Ret g] - 1,0[. 

Now, let us consider the integral 



dt 



3 

2A 



n r(t + ±)r(t + §) 



sin(7rf) 



r(t + i) 



(2.65) 



\(z) = — 

v 1 2m 



c+n+ioa 



dt z * — - — - 

c+n-ioc sm(vrf) 



7T r (t + a) r (t + (3) 



(2.66) 



r(t + i) 

where c = Ret G] — 1,0[ and n is a positive integer so that the path of integration in 
(2.66) lies to the right of the poles of V (t + a) V (t + (3). Then, as n — ► oo, we have, for 



arg(^-)| < 7T and oj = a + (3 — 1, 



l"(*)l 



)• 



(2.67) 



This is the result of a lemma taken from [4] , where cases dealing with more general integrals 
than (2.66) are also considered. One deduces from it that, since (2.65) fullfills the required 
condition, 



IR 



2A 



-Ti— d 



0{e~ n n n+d ). 



(2.68) 



Now if, when |A| is small, we choose n = f& + 6q (where ao > and 6q is bounded), then 



2A 
= O 

= a 



( 


3 




2A 



"FT 







ao 

|A| 



bo+d. 

e 



2oo\ 



O 



-b q 
bo+d 

-ao 



a 



2A 
-bo Q 



2ooy° 



e 1*1 « 



TXT 
o 



O 



l > 
FT 



1 

W\ 



-f%- b °- d a _o 
e W 



2a \ W 

"3"; 



i 



(2.69) 



Minimizing the remainder then implies that ao = § and one finally has 



2|A 



O e 2 M 



(2.70) 
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If one chooses for instance bo such that ^ + = [5^] (the integer part of 2pff)j then 
(2.49) may be written as 

v k=o y 1 

as |A| — > in |arg (^) | < 7r. 

For this optimal value of the perturbative truncation index the remainder is there- 
fore non-perturbative (exponentially suppressed with respect to the coupling constant |A|). 
Eq. (2.71) is called the superasymptotic expansion of Z(0), or optimally truncated hyper- 
asymptotic expansion at zeroth level. 

Now, the remainder (2.65) has been expanded, see (2.60), and a similar analysis may 
be done to deal with the remainder R n , m of this expansion. Indeed, (2.60) reads 



1 m_1 1 
Rn = 7=- (-l)M,- — 



dt 



1m Jd+n-ioo \2XJ sin(Trt) 



7T 



-T(t - j) + R n , m , (2-72) 



with 



1 1 re+n+ioo / 
y/2ir 2i7r J e+n _ ioo \' 



vr 1 



2A / sin(7ra) 2mJ c+m 



c+m+ioo 



dsf(s)T(u- s), (2.73) 



where we recall that f(s) = E(f±lMf±ffi s) . 

m ; r(i)r(|) 

From the same lemma as before and another one that may be found in [4], one finds 
that when n and m become infinite (with c = d and n> m) 

3 



IR, 



2A 
- O 



( 


3 




2A 



e- n (n - m) n - m m m+c - 



(2.74) 



Since we found that n = ^pq- + &o was the optimal value for the truncation index of the 
perturbative series, one may keep this value for n and choose m = ^ + b\ (where a\ > 
and b\ are bounded), so that 



o 



|A|" 



e-la? 1 



3\ 5 /3 \ - 
- a\ 



-ai\ |A| 



as |A| -> in |arg (^) | < it. 



Minimizing the remainder implies a\ = |. This gives 



o 



1 _ 3 1 

p 2|A| 

1 c 3 







2 2 ^ 



|A|-2 



-5^(l+ln2) 



(2.75) 



(2.76) 
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One may therefore write (choosing bo and b\ such that 2m + &o = lw\\] an< ^ STXT + ki = tlpd) 



L 2| A 



Z(0) 



7T 



£ big r (g + 2fc ) A * 



fc=0 



fe! 



(3 



v/2tt 
+ 



JA|- 



3=0 
1 _^L_ 

— e 2 l A l 



(l+ln2) 



(2.77) 



which constitutes the optimally truncated hyper asymptotic expansion of Z(0) at first level, 
in the optimal truncation scheme that we call OTS1 in the following. 

In fact, in OTS1, eq. (2.72) constitutes an asymptotic expansion which, being extracted 
from the exponentially suppressed remainder R 3 , . , allows us to conclude that at the 
first stage of the OTS1 hyper asymptotic approach, we extracted from the perturbative 
remainder a (leading) non-perturbative asymptotic series. 

Let us stress that at first hyper asymptotic level there exists another optimal truncation 
scheme that leads to even better results. Indeed, instead of keeping the superasymptotic 
value n = ^ + bo, one could have let it free: if we choose n = ^ + bo and m = ^ + b\ 
(where oq > 0, a\ > and bo and b\ are bounded), one has 



R n o 



o 



|A|" 



—0.0 



(a - ai) a »- 



■ai 



(2.78) 



as |A| — > in |arg (^) | < ir 

r 

gives 



The minimum of the remainder in this case is obtained for ao = 3 and a\ = |. This 



O 



which is more exponentially suppressed than (2.76). 



(2.79) 



From this, one deduces that (choosing bo and b\ such that ffy +&o = [ ffy ] an d 2pvf +^1 



Z(0) 



tF ^ fc! (3!)* 



v^vr 1 ^ J 2ivrL +[ ^ ] _ loo \3 sin ^ 1 " 



|A|- 



(2.80) 



which constitutes the first level of the hyperasymptotic expansion of Z(0) in the optimal 
truncation scheme that we call OTS2 in the following. 
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Comparing (2.71) and (2.80), one sees that the first level of the hyperasymptotic ex- 
pansion in OTS2 therefore leads to an exponential improvement of the superasymptotic 
expansion. This exponential improvement is of the same order than the exponential im- 
provement that was obtained by the superasymptotic expansion of the perturbative series. 
In this sense eq. (2.80) is non-perturbative. 

The non-perturbative results (2.77) and (2.80) have been obtained from truncation 
indices n = [?&] and m = [j^r] and we saw that the values of the coefficients ao and a\ 
depend on the chosen optimal truncation scheme. 

It is interesting to note that although (2.80) gives a more accurate description of 
Z(0) than (2.77) since the remainder in OTS2 is more exponentially suppressed than the 
remainder in OTS1, the number of perturbative terms retained in (2.80) is at least twice 
the number of those retained in (2.77), which means that taking into account perturbative 
terms that make the perturbative series begin to diverge from the exact result leads to an 
improvement 10 ! 

2.2.4 hyperasymptotics at higher levels and non-perturbative asymptotic re- 
finements 

Our starting point to get the hyperasymptotic expansion of Z(0) at second level is the 
battle-horse formula (2.40) obtained by the formal strategy developed in section 2.2.1 and 
which has been rigorously demonstrated in section 2.2.2, namely 

n-i , „ x -k 



Z (o) = B-iM(s) 

fc=0 v 7 

J ^Jd+n-ioc V2A; sin(Trf) V J> 

+ Rn,m , (2-81) 



where 



i i re+n+ico / o \ — u i rc+m+ioo 

R n , m = --=- — / dul — ) ^~r^w- / dsf(s)T(u - s), (2.82) 

V27T ^*7T Je+n-ioc \ ZA J Sm{lTU) ZITT 

with n > m, c e] — 1,0[ and d g] — 1,0[. In the double Mellin-Barnes integral, e G 
] — min(l, n — m — c), 0[. 

The second hyperasymptotic level will appear from an expansion of the double Mellin- 
Barnes remainder integral R n ,m- We do not give the details of the calculations since they 
are performed in exactly the same manner as what we did on the Mellin-Barnes remainder 
integral (2.57) in section 2.2.2. 

Noticing that 

r(5+J)l> + i) vr 



m = - 



r(i)r(f)r(s+l) sin(vr S ) 



3 This will be confirmed numerically in section 2.2.5. 
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TT 



1 



sin(7r S ) r (I) r (|) 



h+m' +ioo 



m'-l "~ r "' 

^(-i)%r( s -/)+ | 



eft 
2ivr 



/(t)r( s -t) 



h+m' —ioo 

where /i + ml < Res and h G]0, 1[, the final result reads 



(2.83) 



n— 1 



fc=0 



Z(0)=£(-1)M* 

_L Y(-iyA / 



d+n+ioo 



i=o 



d+n-ioc V 2A / sin(TTt) 



:r(i-i) 



m — 1 



2ivr 



2 re+n+ioo 



du 



e+n—ioo 



TT 



2A / sin(7r«) 

c+m+ioo 



+ 



(V27T) 



rc+rn+ioo 

x / dsT(s-l)— — -r(n-s) 

Jc+m-ioo Sm(7TSj 
: / 1 \ 3 f e + n +i°° / 3 \ ~ W TT rc+m+ioo ^ 

z — 9 ( ttt- I / du I — I -— — - / ds -— — rT(n — s) 

2tt) V 2 W Je+n-ioo V 2A / sm (™) Jc+m-ioo Sin(TTS) 

>-/i+m'+ioo 

x / tft f(t)T(s-t). (2.84) 

J h+m' —ioo 



The constraint /i + m' < Re s implies /i + ml < c + m in (2.84). 

We have now three interwoven partial sums where the third one is expressed in terms 
of so-called hyperterminants [4] (here defined as double Mellin-Barnes integrals). 

For the iV-point function case, we have, at second hyper asymptotic level 11 , 



,2p -i0 2 -A^4 



2 2p (\ p\ _V3 p 



/+oo 
-oo 



\k D (-P) 



2A 



-n — p + j + 1, 



2A 



\i r(-p) 



c+m+ioo 



J 2m \2\) sin(Trs) 



TT 



c+m— ioo 



r(s — p — l)T(n + p — s) T ( — n — p + s + 1, 



2A 



_ ( AY 



e 2A 



c+m+ioo 



r(i-I)r(|-I) 

3 x - s 



TT 



c+m— ioo 



ds 



2m \2X J sin(7rs) 



r(n + p — s) T — n — p + s + 1 



2A 



11 We give here the hyperasymptotic expansion in an alternative form where (hyper)terminants are written 
in terms of (integrals of) incomplete gamma functions. 
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c+m +ioo 

X 

c+m'—ioo 

with 



/dt 
— g_ p (t)T(s-p-t) (2.85) 



9p(s) _ r(. + i- S )r(. + i- S 



TT 



r (3 - 5) r (! - 1) r( s + 1) sinM 

7T 1 

S in(^)r(i-f)r(f-|) 

" m '_l c+ra'+ioo 



2 (-i)^;- p) r( s - p - + y — r( a - v - 1) 

1=0 c+m'-ioo 



(2.86) 



The resurgence phenomenon is clearly apparent in (2.84) since the coefficient Ak ap- 
pears at each hyper asymptotic order. One may however notice that for the iV-point func- 

(v) (—v) 

tions case, resurgence manifests itself in a different way: the coefficients Bj ' and B K - 
appear alternately at successive levels in the expansion. 

This difference comes from the fact that for Z(0) the function f(s) in the l.h.s. of 
(2.83) appears also under the integral sign in the r.h.s., whereas for the iV-point functions 
the function g p (s) in the l.h.s. of (2.86) appears as g~ p (s) under the integral sign in the 
r.h.s. 

Thanks to these interesting symmetries that are reflected by the resurgence phe- 
nomenon, one may deduce that in our cases of study, hyper asymptotic expansions at 
arbitrary hyperasymptotic levels may be obtained in a straightforward way. 

Let us now discuss optimal truncation schemes at second hyperasymptotic level. 
In OTS1, one finds that the optimal values for the truncation indices are n = [gpq-], 
[ifXf] an d m' = [gpq-L while in OTS2 these values are n = [^y], m = [pq-] and 



m 



2.2.5 Numerical analysis 

In this section we want to compare numerical results that may be obtained from the 
hyperasymptotic expansion of Z(0) at various hyperasymptotic levels, in both cases of 
optimal truncation schemes OTS1 and OTS2, with those obtained from perturbation theory 
in section 2.1.3, i.e. for A = |. 

Recall that from (2.7) one gets, with an 8 decimal places precision, the "exact" result 

Z(0)| A= i = 0.96556048, (2.87) 

while the best perturbative prediction is given by (see (2.19)) 

Z(0)|^*i = 0.96555187 ± 0.00140990. (2.88) 



First hyperasymptotic level 
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From (2.40) (or (2.34)) in OTS1 without considering the remainder integral 12 , i.e. 
from (2.77), one obtains 

^(O)!^ 3 ! 1 = 0.96552297, (2.89) 

A— 3 

while in OTS2, i.e from (2.80), one has 

Z(0)|°ZT = 0.96556291, (2.90) 

3 

Comparing with (2.88) one sees that perturbation theory gives a better prediction than 
(2.89). This is due to the fact that the perturbative contribution to (2.89) comes from a 
truncation of (2.8) after the fourth term, while (2.88) is not exactly a naive truncation of 
the perturbative series (see (2.21)). (2.90) is however closer to the exact result than the 
central value of (2.88). 

Second hyperasymptotic level 

From (2.84) in OTS1 (i.e. with n = [^y], m = [^] and rri = [g^]) still without 
considering the remainder integral, one gets 

Z(0)\?if = 0.96556492, (2.91) 

A— 3 

while in OTS2, i.e with n = [^], m = [^] and m! = [^], 

Z(0)\°l¥ = 0.96556046, (2.92) 

A— 3 

At this hyperasymptotic level (2.92) is already very close to (2.87) while (2.91) gives a 
prediction less good than (2.90) but better than perturbation theory. 

Third hyperasymptotic level 

At this level (for which we do not give the expression of the corresponding hyper- 
asymptotic expansion of Z(0)), one has no prediction from OTS1, since it would imply 
m" = [jgkj] = (where m" is the truncation index of the partial sum appearing at third 
hyperasymptotic level) . 

In OTS2, however, one has n = [r|r], m = [^pq-], m' = [j^] and m" = [jpq-], so that 

Z(0)|°3 2 = 0.965560486, (2.93) 

A— 3 

which gives the exact result with an 8 decimal places precision, although we see that n = 18, 
which means that the perturbative partial sum contribution to (2.93) is (see Table 1) 

= -27.6964871, (2.94) 

a value extremely far from the exact result (2.87)! 

It is clear from this analysis that OTS2 is a much better optimal truncation scheme 
than OTS1 and that the results obtained in OTS2 give much closer value to the exact 
result than what perturbation theory does. 

12 In all the numerical evaluations coming from hyperasymptotic formulas, adding the remainder integral 
of course always leads to the exact result (2.87). 
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3. Conclusions 



This paper was aimed to show how an interpretative procedure of divergent perturbative 
series, involving Mellin-Barnes representation and Borel resummation, leads to a non- 
perturbative asymptotic improvement of the iV-point functions perturbative expansions, 
for the simple case of zero-dimensional euclidean 4 theory, by the appearance of hyper- 
asymptotic expansions. We saw that these hyperasymptotic expansions are composed of 
interwoven partial sums whose coefficients, in our cases of study, are linked together by a 
resurgence phenomenon. The non-perturbative interpretation of our results relies crucially 
on so-called optimal truncation schemes of these partial sums. A numerical analysis has 
been performed, showing the much more accurate results that one may obtain with the 
non-perturbative hyperasymptotic expansions, compared to the traditional perturbative 
approach. We also saw that on the Stokes line A < 0, one may find, from hyperasymptotic 
expansions, imaginary contributions that are not obtainable from the truncated perturba- 
tive expansions (which are non-Borel summable for these particular values of the coupling 
constant). 

One of the important conclusions concerns the two different optimal truncation schemes 
that were studied in this paper. Indeed, the best one (i.e. the one that leads to the best 
predictions) implies a truncation of the perturbative series that depends on the level at 
which the hyperasymptotic improvement is performed. This leads to the striking result 
that the corresponding (best) predictions, which are of an amazing accuracy (see eq. (2.93) 
compared to the exact result (2.87)), involve perturbative contributions that, taken inde- 
pendently of the non-perturbative corrections, would lead to desastrous results (see eq. 
(2.94)). This is due to the fact that the more one wants to increase the exponential 
improvement (i.e. the more one wants to include non-perturbative contributions in the 
analysis), the more one also has to include perturbative contributions far in the divergent 
tail of the perturbative series. 

These results clearly show that hyperasymptotic expansions are tools that could have 
a lot of applications in high energy physics. In this respect, concerning the study of 
Stokes phenomenon (and the smoothing of the Borel ambiguity) , we would like to mention 
an interesting paper [12] which appeared a few days before the end of the writting of 
our manuscript and dealing, among others things, with a hyperasymptotic approach of 
instantons in topological string theory and c = 1 matrix models. It is however not based 
on Mellin-Barnes hyperasymptotic theory but rather on the hyperasymptotic approach 
developped by Berry and Howls [13]. 

Acknowledgements We would like to thank Santi Peris and Eduardo de Rafael for their 
comments on the manuscript. D.G. acknowledges financial support from CICYT-FEDER- 
FPA2008-01430, and the Spanish Consolider-Ingenio 2010 Program CPAN (CSD2007- 
00042). 
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A. General bound for the remainder of an alternating asymptotic series 

Let = J2T=o Uk ^ e the sum °f a divergent alternating asymptotic series. In the inter- 
pretative procedure to give a meaning to this sum, one defines the remainder R n as 

5*oo = S n -i + R n . (A.l) 

We therefore have 

Rn — Rn+i = u n , (A. 2) 

so that if u n > 0, since the asymptotic series is an alternating one, R n +i < and then 
Rn > 0. Therefeore R n - R n +i = \R n \ + \R n +i\ = \u n \. 

On the other hand, if u n < 0, then R n +\ > and R n < 0. Therefeore R n — R n +i = 

— \R n \ — |.Rn+l| = — \u n \- 

We conclude that in all cases \R n \ + |-R n +i| = \u n \ and that we therefore also have 
|-Rn-i| + \Rn\ = \u n -i\, which means that 

\R n \ < \u n \ (A.3) 

and 

|-Rn| < |«n-l|- (A.4) 



B. Inverse factorial expansion and Barnes Lemma 

In this appendix, we show how to prove the inverse factorial expansion eq.(2.25). We also 
show how to find (2.57) from (2.33). 

Barnes lemma, applied to our case, reads [4] 

r (k + \) r (k + 1) _ i r c +™ r ( s + \) r ( s + f ) r(k - 8 )r(- 8 ) 



i p r ( S + j) r [s + 1 r(fc - 8 )T(- 8 ) 
-Loo ds « ' 



r(fc + i) " 2ivr y c _ 1QO r(i)r(|) 

where c = Res e] — \ ,0[ and A; is a positive integer with k > n where n is the truncation 
index of the perturbative series. 

By Stirling formula (2.50) and Cauchy's theorem, one may prove that 

T(k + \)Y(k + \) _ Resid+ ± p^ A r( J + i)r( J + §)r( fc - a )r(- J ) > 



r(fc + i) " ' 2ivr J d+m _ loo rQ)r(|) 

where Resid. denotes the sum of the residues of the poles of the integrand, located in the 
interval ]0, d + m[ (with d €.] - 1, 0[). 

Looking at the integrand, it is clear that in the considered interval, only T(k — s) and 
r(— s) may have poles. Here comes the crucial point. We do not want that T(k — s) gives 
contributions. For this we have to impose that m < n, since k > n. With this constraint 
one gets only the contribution of the residues of r(— s) and this implies the inverse factorial 
expansion 

1 rf! -E^W-il+r / dsf(snk-s), (B.3) 
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where 

rhy r (s + 1) r („ + 1) r(-») 

m mm — ( ' 

and 

One may prove that (B.3) is an asymptotic expansion when k (or n) goes to infinity [6]. 

If one does not impose the constraint, i.e. if m > n, then one has the alternative 
expansion 

r n^i) +l) -B-^-,) 

V 1 3=0 



m-1 , r /l \ /3 



A J T 14 



+ — / /( S )r(fc - s) , (B.6) 

which cannot be at the basis of the hyperasymptotic procedure that we show in this paper. 



As we said at the beginning of this appendix, we will now demonstrate how to find 
(2.59) from (2.33). 

Recall that (2.33) reads 

1 ^ (-l) fc r(l + 2fc) 
k\ (3!) fe 

k=n 

-^g(-W-^)(|)"' r (-" +J ' + 4) 

( — ~\\ n s 1 rc+m+ioo /o\-s / o \ 

+ ht c ^L^ d ' («) /(s,r(n - s,r (- + s + 11 5a) • (R7 > 

where c e] — 1, 0[. First we choose m = since (2.59) is independant of m. This gives 

i ^ (-i)*r(^2fc) k 

fc! (3!) fe 



(-l) n JL 1 
-e2A . 



pc+ioo / o \ -s / o \ 

L/'U) /Wr(»-.)r(-» + . + i, s ), P» 



v^vr 2i7T 

where now c e] — |,0[ (see the comment after (2.25)). 

Then we need the Mellin-barnes representation of the incomplete gamma function, eq. 
(2.61), that we demonstrate now. 

(o \ poo / o \ s— n poo / o\ \ s— n 

- n + s + l,-)=l dy e- = e"« (^) jf ^ (l + y *) e-. 

(B.9) 
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Using 

d+ioo 

J 2m Y(n — s) 

d—ico 

which is valid for d = Ret g]0, n — Res[ (where Res = c) and | argX| < |, one has 

d-\-ioo , 

3\_„-AfZ\ S ~ n 1 f ,j2\y t T(t)T(n-s-t) r^„-t„-x 



r -n + s + 1, — =e~5X — / di — v - , -/ dx x~ l e 

'2XJ \2\J 2m J V 3 J r(n-s) J 



d—ioo 



d+ioo 

1 s / 3 V~ n 1 /" , /2A\~* vr 



(B.H) 

d— ioo 

(B.12) 

where the y-integral imposes Ret < 1, so that since n > and c < 0, d is now restricted 
to belong to ]0, 1[. 
Therefore, 

i ^(-i) fc r(i + 2fc) 

fc=n 

c+ioo (i+ioo 

(-l) n /3 V" f ds „ , /" dt /2A\~* vr . 

' ' ' /(*) / 77- hr -^-rr(n- S -t) (B.14) 



v^vr \ 2 V 7 2i7r J 2iir \ 3 J sm(irt) 

c—ioo d—ioo 

Since c g] — |,0[, by Barnes Lemma, one has 

c—ioo 

Inserting this expression in (B.14), one obtains 

i ^(-i) fc r(i + 2fc) 

V^nti (3^ A ( 6) 

k=n 

d+ioo 

_(-!)" /3\" /• dt /2A\- f tt rQ + n-t)r(f + n-t) 

^2vr V 2A / 7 2ivr V 3 / sinvrt r(l + n - 1) { ' ' 

d—ioo 

Performing the change of variable s = n — t, one finds 

i +g(-i)fc r(i + 2AQ fc 

V^e^i (3l)^ A ( 8) 



fc=n 

c+n+ico 



/ fir 8 i r (^+i) r (^+i) fB191 

v^vr 7 2m \2\) sin^s 1> + 1) K ' ' 



c+n—ioo 
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where c g] — 1, 0[. This last equation is (2.59). 

It is clear that one could also find this result by starting the resummation of the tail 
(2.24) with Barnes lemma instead of the inverse factorial expansion (2.25) (since we see 
that putting m = in (2.25), which is nothing but Barnes lemma (B.l), is equivalent to 
have m = in (2.33)). 
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